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ABSTRACT: We extend our study of all-order linearly resummed hydrodynamics in a flat 
space [1, 2| to fluids in weakly curved spaces. The underlying microscopic theory is a 
finite temperature J\T = 4 super-Yang-Mills theory at strong coupling. The AdS/CFT 
correspondence relates black brane solutions of the Einstein gravity in asymptotically locally 
AdSs geometry to relativistic conformal fluids in a weakly curved 4D background. To 
linear order in the amplitude of hydrodynamic variables and metric perturbations, the 
fluid’s energy-momentum tensor is computed with derivatives of both the fluid velocity 
and background metric resummed to all orders. We extensively discuss the meaning of all 
order hydrodynamics by expressing it in terms of the memory function formalism, which 
is also suitable for practical simulations. In addition to two viscosity functions discussed 
at length in refs. [1, 2], we find four curvature induced structures coupled to the fluid via 
new transport coefficient functions. In ref. [3], the latter were referred to as gravitational 
susceptibilities of the fluid. We analytically compute these coefficients in the hydrodynamic 
limit, and then numerically up to large values of momenta. 
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1 Introduction 

Fluid dynamics [4-6] is an effective long-distance description valid for most classical or 
quantum many-body systems at nonzero temperature. It is defined in terms of thermody¬ 
namical variables. Derivative expansion in fluid-dynamic variables such as velocity accounts 
for deviations from thermal equilibrium. At each order, the derivative expansion is fixed by 
thermodynamic considerations and symmetries, up to a finite number of transport coeffi¬ 
cients, such as viscosity and conductivity. The latter are not calculable from hydrodynamics 
itself, but have to be determined from underlying microscopic theory or experimentally. 

The stress-energy tensor T^ u of a relativistic fluid is conveniently written as 

Tf_u; = (e T P'jUfxUi/ + Pg^u T (1-1) 
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where e, u^ and g^ v are the fluid energy density, four-velocity and metric tensor. The 
microscopic theory underlying the hydrodynamics is supposed to specify the fluid pressure 
P via equation of state, P = P(e). The viscous effects are encoded in the dissipation tensor 

IV> EE {u a/3 + u Pa ) - ^A M „A^n a/3 , (1.2) 

where, in the local rest frame, = g^ + u^u^ acts as a projector on the spatial subspace. 
For conformal fluids, the first order gradient expansion has one term only, the Navier-Stokes 
term 

= —2'^oV^tty, (1-3) 

where ?yo is a shear viscosity. Eq. (1.3) relates the spatial components of the stress-energy 
tensor (Tjj) to time involving components T$j. A relation of this type is frequently referred 
to as a constitutive relation. Energy/momentum conservation (conservation of T^ u ) leads 
to fluid dynamical equations, which are normally solved as initial value problem, when 
initial profiles for temperature and velocity are specified. We shall generically refer to these 
equations as Navier-Stokes equations, even when we extend the discussion beyond the first 
gradient. 

Although fluid dynamics has long history, theoretical foundations of relativistic viscous 
hydrodynamics are not yet fully established. The Navier-Stokes hydrodynamics introduced 
in (1.3) leads to violations of causality: the set of fluid dynamical equations makes it 
possible to propagate signals faster than light. To overcome this problem, simulations of 
relativistic hydrodynamics are usually based on phenomenological prescriptions of [7-10], 
which admix viscous effects from second order derivatives, so to make the fluid dynamical 
equations causal. Refs. [7-10] introduced retardation effects for irreversible currents, which, 
via equations of motion, become additional degrees of freedom. In other words, one needs to 
include higher order gradient terms in the derivative expansion in order to obtain a causal 
formulation. In general, causality is violated if the derivative expansion is truncated at any 
fixed, order. It is supposed to be restored when all order gradient terms are included, which 
we refer to as all order resummed hydrodynamics. Causality usually also implies stability 
and even can be used to constrain possible values of higher order transport coefficients [11— 
16]. 

Discovery of AdS/CFT correspondence [17-19], and its generalizations, has established 
a novel approach for exploring a large class of quantum field theories at strong coupling. 
Within the holographic picture, hydrodynamic fluctuations can be regarded as gravitational 
fluctuations of black holes in asymptotic AdS spacetime, and vice versa [20, 21]. The 
connection between fluid dynamics and black holes in AdS gravity is usually referred to as 
fluid/gravity correspondence. One central result of this correspondence is the ratio between 
r/o and the entropy density s [20, 22, 23] 


rjo _ 
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This ratio (1.4) is universal [24-26] for a large class of strongly interacting gauge theories for 
which holographic duals are governed by Einstein gravity in asymptotically AdS spacetime. 


Remarkably, the fluid/gravity correspondence is not limited to linear response theory 
for small perturbations of the velocity held. Nonlinearly generalized Navier-Stokes equa¬ 
tions are completely encoded in the Einstein equations. Particularly, the formalism of [27] 
provides a systematic framework to construct nonlinear fluid dynamics, order by order 
in the velocity derivative expansion, with the transport coefficients determined from the 
gravity side. The study of [27] was subsequently generalized to conformal fluids in higher 
dimensions [28], weakly curved background manifolds [29], and to forced fluids [30]. We 
refer the reader to [31-33] for comprehensive reviews of fluid/gravity correspondence. 

In [1, 2], we built upon previous work of [3] and constructed a hat space all-order linearly 
resummed relativistic conformal hydrodynamics using the huid/gravity correspondence. We 
have collected all the derivative terms that are linear in the fluid dynamic variables, like 
(VV • • • 'S7u) liu ] but we neglected all nonlinear structures, such as (Vu)^. The velocity held 
u^ has a big time component, uq ~ 1, whereas the three velocity u is considered to be small. 
It is worth stressing that the linear derivative terms can be made dominant by restricting 
our study to small u, hence our approximation is well under control. Ideologically similar to 
our study here, resummations of gradient terms for boost-invariant plasma were considered 
in [34-39], with the prime difference that non-linear terms were also accounted for in these 
papers. 

In [1, 2] relativistic hydrodynamics with all order derivatives resummed was found to 
have a rich structure, absent in a strict low frequency/momentum approximation. The 
fluid stress-energy tensor was expressed using the shear term (1.3), with 7]q replaced by 
a viscosity function of space-time derivative operators, and a new viscous term, which 
emerged starting from the third order in the gradient expansion. In Fourier space, both 
coefficient functions accompanying these two terms become functions of frequency and 
spatial momentum. They were calculated in [1, 2], Those viscosity functions were found to 
vanish at very large momenta, which is a marking signature of causality restoration. 

In the present work, we generalize computations of [1, 2] by including small metric per¬ 
turbations in the boundary theory. For A/" = 4 super-Yang-Mills theory, we deduce the stress 
tensor for all-order linearized relativistic hydrodynamics in a weakly curved non-dynamical 
background spacetime. Following the method of [27, 30], we construct asymptotically locally 
AdSs solutions to the Einstein gravity with a weakly curved boundary metric. However, 
our procedure is somewhat different from that of [27, 30], and it has been invented in order 
to include all order linear structures in a self-consistent manner. In particular, rather than 
considering order-by-order derivative expansion, we will collect all the derivative terms in a 
unified way [1, 2], We first solve the dynamical components of the bulk Einstein equations, 
which is sufficient to compute an “off-shell” stress-energy tensor of the boundary fluid with 
the transport coefficients to be determined completely. Putting “on-shell” thus obtained 
stress-energy tensor is equivalent to constraints of the Einstein equations. 

The boundary metric perturbation can be regarded as an external force acting on the 
fluid in flat space [30]. Consider the covariant Navier-Stokes equations 

= 0, (1.5) 

where is compatible with g /iL/ . For a weakly curved spacetime, g /lL/ = rj^ + h^ u , where 
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r)^ is a Minkowski metric. For the rest of this paper we keep terms up to the first order in 
hfj, u only. The Navier-Stokes equation (1.5) can be rewritten as 


d^To)^ = ( 1 . 6 ) 

where (To) /ti , is stress-energy tensor of the fluid in Minkowski space. The effective forcing 
term f v is a functional of and Uu. Its explicit form can be found in [30]. By appropriately 
choosing h )w one can stir the fluid into various flows. 

We now present the main results of this paper, before embarking on their derivation. 
The fluid energy density is related to pressure via equation of state 

e = 3P = 3(7tT) 4 , (1.7) 

consistent with the tracelessness condition T/, 1 = 0. The tensor II^j, is 

If hi/ — T ku T pu Cjiav ft . . 

+ ^V a V^C m ^-6u a V a R 

where C liau p and R lw are the Weyl and Ricci tensors of g^. Appendix A (eq. (A. 4)) 
provides expressions for these tensors in terms of the metric derivatives, assembled into 
useful combinations listed in Table 1. The transport coefficients are functions of derivative 
operators 


r/Kv ai v a v a ], CKv a ,v a v a ], «Kv a ,v a v“], 

(Kv a ,v a v“], 0Kv a ,v a v a ]. 

In Fourier space, via the replacement d^ —> (—iu,qi), these operators become functions of 
u ; and q 2 . Below, we will be frequently using mixed notations: the tensor structures 
etc. are usually written explicitly in the form of derivatives while the transport coefficients 
will be be always presented as functions of momenta. The constitutive relation (1.8) was 
introduced in [3], but with £- and 0-terms missed. 

The reader might be puzzled by the fact that the constitutive relation (1.8) contains 
infinitely many time derivatives. As a result, the all-order resummed hydrodynamics at 
hand does not seem to be solvable as initial value problem, because infinitely many initial 
conditions would need to be supplied. To some extent this is indeed true. Nevertheless we 
can make sense of it by transforming the gradient expansion of the transport coefficients 
into memory functions with a well-prescribed initial value setup. In section 2, we illustrate 
our idea on a simple model of diffusion in magnetization. 

The viscosity functions r/ and ( were computed in [1, 2], and the main results will be 
quoted below. The prime focus of the present paper will be on the rest of the transport 
coefficients that couple the stress tensor to the curvatures. In [3] these coefficients were 
named gravitational susceptibility of the fluid (GSFs). GSFs measure fluid’s response to 
external gravitational perturbations. In the hydrodynamic limit, uj —»• 0 and q —> 0, the 
viscosities and the GSFs can be expanded in powers of momenta, and we were able to 
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compute a few terms in the expansions analytically 

r)= 1 + ^(2 —ln2 )iu- ^g 2 - ^ [6vr - vr 2 + 12 (2 - 31n2 + In 2 2)] w 2 + • • • , 

C = ~ (5 — VT - 2 In 2) H-, (1.9) 

k = 2 + ^ (5 + 7T - 6 log 2) iu 4-, p = 2 4-, £ = In 2 — ^ 4- , 6 = 

The results are presented in terms of dimensionless frequency uj and spatial momenta gj. For 
the rest of the paper the choice of units is set by nT = 1. The physical momenta should 
be understood as nTu and irTqi. The first term in n was computed in [29, 40]. p was 
introduced in [3] and agrees with the result here. The relation 6 = 3£/2 holds for generic 
values of momenta. We complete the study of the GSFs by evaluating them numerically 
up to very large values of momenta and report relevant results in subsection 4.3. We also 
present a time dependent memory function derived from the shear viscosity function r]. For 
convenience of possible applications, all our numerical results for the transport coefficient 
functions are deposited into a Mathematica file downloadable from the link [41]. 

The technical parts of this paper are presented in sections 3 and 4, which follow the 
setup of [1, 2], To facilitate reading of these parts, we now briefly sketch the main steps of 
our calculation. In section 3, we deal with the gravitational dynamics of a weakly perturbed 
boosted black brane in AdSs- The perturbation is set to be driven by varying boundary 
velocity held u^{x) and temperature T{x). The boundary conditions on the 5D bulk metric 
are imposed so that the boundary metric is rj^ + h^. We then parameterize the bulk 
metric perturbation in terms of ten functions h, k, ji and cx{j . which are both functions 
of the holographic coordinate and functionals of u^ and h^ u . We then holographically 
read off the boundary stress tensor and express it in terms of, at this stage yet unknown, 
near-boundary asymptotic behavior of h, k, ji and a t j. 

In section 4, we solve the Einstein equations for h, k, ji and aij. This is done in two 
steps. First, from the vector and tensor h IJV we construct all possible linear structures 
(scalars, vectors, and tensors), which are summarized in Table 1. These structures serve as 
a basis for linear decomposition of h, k , ji and aij. The coefficients are generally denoted as 
Vi and Tj. which depend on four-momenta and the holographic coordinate only. What we 
achieve is that the functional dependence on unknown velocity field and boundary metric 
is completely removed, and we are left with ordinary second order differential equations for 
Vj and T, l . Second, we solve these equations. The results for the viscosities and the GSFs 
are obtained from a pre-asymptotic, near the boundary behavior of V) and Tj. Summary 
and discussion can be found in section 5. We deposited many computational details into 
several Appendices. 

2 All order resummed hydrodynamics 

In this section we provide a clarification about what we actually mean by all-order resummed 
hydrodynamics. As has been mentioned in the Introduction, in the constitutive relation 
for the dissipation tensor we resunnn all orders in the gradient expansion, including infinite 
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number of time derivatives. This particularly means the effective dynamical equations in 
principle require an infinite set of initial conditions or, equivalently, the hydrodynamics at 
hand is a theory of infinite number of degrees of freedom. Indeed, these are the quasi-normal 
modes of the dual gravitational theory. An alternative way to see the origin of the problem 
is to realize that exact gravitational theory in the bulk cannot be mapped onto a single 
degree of freedom on the boundary. Below we are to solve for the bulk perturbations as a 
boundary value problem, with boundary conditions imposed both at the black hole horizon 
and the AdS boundary. If, to the contrary, we were to solve the bulk dynamics as initial 
value gravitational problem, we would have to supply initial conditions for the gravitational 
perturbation, but in entire bulk. Roughly speaking, the initial condition along the extra 
dimension (holographic coordinate) would have to be mapped into infinitely many initial 
conditions on the boundary. 

Yet, it turns out that it is possible to formulate this all-order hydrodynamics as a 
normal initial value problem with all the time derivatives absorbed into a complicated 
memory function. We illustrate our formalism by focusing on a classic problem of spin 
magnetization in external field [5, 6], which is free of unnecessary complexities related to 
relativistic formalism and tensorial structures, but otherwise conceptually identical to the 
hydrodynamics we would like to study. 

Consider spin magnetization (in some direction) M(t,x ), which can be created by an 
external magnetic field H(t, x ). We also define the magnetization current J(t, x). Compared 
with the hydrodynamics, M is identified with the velocity H plays the same role as 
the external gravitational held, and J is analogous to 11^. Similarly to our hydrodynamic 
construction, for the current we introduce the constitutive relation 

J(t, x ) = £>(<%, V 2 )VM(t, x) + a(d t , V 2 ) H(t, x). (2.1) 

Just like the viscosities and GSFs, the diffusion coefficient D and the magnetic susceptibility 
a are considered to be functionals of time and space derivatives. In Fourier space this reads 

Jiiu^q) = D(uj,q 2 )qiM(uj,q 2 ) + q 2 ) H(uj, q 2 ). (2.2) 


The continuity equation 

M + V • J = 0 (2.3) 


defines the dynamics of the system just like the conservation law for the energy momentum 
tensor while the constitutive relation (2.2) is analogous to (1.8). Fourier transforming (2.2) 
back, we can have a different representation of (2.1) (suppressing the spatial coordinates): 



dt' 


D(t 


f')VM(f') + ef(t - t') H(t') 


(2.4) 


Here 

D(t — t') = [ — D{uj) e -^-^ ■ cr(t-t') = [ — a ( u ) (2.5) 

J 27t J 2-1 r 

The current J at time t should be affected only by the state of the system in the past. This 
causality requirement implies that both response functions have no support in the future: 
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D(t ) ~ 0(t), and a(t) ~ 0(i). This is achieved only when the poles of D(u>) and <t(w) are 
all lying below the real axis. This must be the case for any causal theory. Thus 

rt 


J(t) = 


d£ 


D(t 


f')VM(f') + er(t - t') H(t') 


( 2 . 6 ) 


The external magnetic field H is normally turned on at negative times, so to create 
initial magnetization at t = 0, and then turned off at t = 0 ( H(t > 0) = 0), letting the 
system to freely relax to its equilibrium at infinite future. For such experimental setup, for 
positive times the current 


J(t > 0) = / dt'D{t - t')VM(t') + J H (t) 


(2.7) 


with 




d£ 


D(t - t')VM(t') + a(t - t') H(t') 


( 2 . 8 ) 


Generically, Jh is not vanishing and it accounts for the entire history of the system at 
negative times. This is in contrast to a typical memory function-based approach, where 
one introduces constitutive relation (2.7) assuming Jh = 0 and then also models D [5, 6]. 

Our construction is so far formally exact. However, in order to solve the dynamical 
equation (2.3), it is not sufficient to provide the initial condition for magnetization only, 
but we also need the “history” current Jh at all times, equivalent to providing infinitely 
many additional initial conditions. 

We are now to discuss under what conditions we can nevertheless set Jh to zero, 
casting our theory into a well-defined initial value problem. The response functions D and 
a are some given functions defined by underlying microscopic theory. Thus the equation 
J#(f > 0) = 0 is in fact an equation for the magnetic field H at t < 0. It is, however, not 
obvious that there exists a solution for generic D and a, because we want the current Jh = 0 
at all positive times. Even though we cannot guarantee vanishing of the current identically, 
it is safe to assume that its effect could be rendered negligibly small. Particularly, it is 
obvious that Jh vanishes at late times and its only potential influence could be at very 
early times when the current J ~ Jh- 

When it comes to modeling of D , a first try is usually the relaxation time approxima¬ 
tion [7-10]: 

D IS (t) = Doe-^K (2.9) 


This model assumes D(u) has a purely imaginary and momentum independent pole at 
uj = —i/tr. A further improvement of this model would be to account for additional poles, 
say, u± = ±wr — iujj [42] : 1 


D±(t) = doe t ^ TR + d\e tuJl cos (ujRt) 


( 2 . 10 ) 


In contrast to the above models, a calculation analogous to one presented below for the 
viscosity functions and the GSFs, would determine D and a exactly, from the underlying 

1 With an equivalent model for the viscosity function, we almost recover eqs. (12),(15),(16) of [42], 
However, our version of eq. (12) has an extra term and a source term. 
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microscopic theory. Any practical application of thus resumed hydrodynamics would have 
to additionally impose Jh = 0 . 

After this explanatory section, we return to the main bulk of our work, and that is to 
determine the viscosity functions and the GSFs for AT = 4 CFT precisely. 


3 Linearized fluid/gravity correspondence: weakly curved boundary 


Notations and Conventions. Upper case Latin indices {M, N, ■ ■ ■ } and lower case 
Greek indices {n, u, ■ ■ ■ } denote the bulk and boundary directions, respectively. Lower 
case Latin indices {i, j, • • ■ } are used to specify spatial directions on the boundary. Indices 
v and 0 will be used alternatively to denote the time direction. We take the AdSs radius to 
be unity which is equivalent to setting the bulk cosmological constant A = —6. The bulk 
metric is Gmn, and denotes (non-dynamical) metric in the boundary theory. Through¬ 
out this paper, we work with the mostly plus signature for the metric tensors. Further 
details can be found in Appendix A. 

A universal subsector of the AdS/CFT correspondence is the Einstein-Hilbert term 
plus a negative cosmological constant, 

s eh = y^^ f d 5 X V^G( R +l2). (3.1) 

Its variation leads to Einstein equations in the bulk, 

Emn = Rmn ~ GmnR ~ §Gmn = 0 . (3.2) 

For the variational principle to be well-defined, the action (3.1) has to be supplemented 
with the Gibbons-Hawking surface term 

Sgr = g7r ^ J d^xy/^y K[ 7 ] (3.3) 

where 7 ^ is an induced metric on a fixed r hypersurface £ with r being the holographically 
emergent coordinate. The extrinsic curvature K^ u of £ is 


E-fii/ — 7 


(3.4) 


where n^ is an outgoing vector, normal to £ and V is compatible with The boundary 
metric g, JiU is related to 7 ^ through g^ = lirn^oo ( 7 ^/r 2 ). 

To ensure finiteness of the boundary stress-energy tensor, we use holographic renormal¬ 
isation and add to the actions (3.1,3. 3) a suitable counter-term S c .t.. For curved boundary 
relevant to our study, the counter-term was worked out in [43, 44] 


s c.t. = ~ iq Iq n J d 4 xy/^y | 6 + ^^[ 7 ] - logl Q«H7]^/47] - ^[7])}, (3-5) 

which is a functional of the induced metric The boundary stress-energy tensor T^ u is 
defined as follows 


T'nv = hm T liv (r) = - lim 

r—¥ 00 1 —>-oo 



2 5S ven \ 
yj —7 J ’ 
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(3.6) 






where the renormalized action S ren is 


^ren = -S'eH + ScH + <5c.t. ■ (3.7) 

From the definition (3.6), we have 

Tfi U (r) = —2 r ~ K'ltiv + 37^ — \ + 2 T fW log —^. (3-8) 

T^(r) T ^R) 

where G^ u {r) is the Einstein tensor compatible with 7 ^. The expression for T“„ can be 
found in [43 ] 2 * 

TfZv = — ^Rfipv\[g\R p [g] + t^V^Vi ,R[g\ — -V R^[g] + —R[g\Rp,u[g\ 

1 (3-9) 

+ 48 g,,u (v 2 R[g] - R 2 [g] + 3 R p x[g]R pX [g}) ■ 

When the boundary is flat, g^ u —>• g^, T° w vanishes and we are left with T )”,. which was first 
worked out in [45]. The second piece, T“ y , is crucial in removing logarithmic divergences, 
which emerge in T™ u due to boundary curvature. Finally, in (3.8) we normalized the Newton 
constant Gn by IGitGn = 1 so that e and P take the forms (1.7). 

In order to incorporate metric perturbation in the boundary theory, we follow [29, 30] 
and consider the boosted asymptotically locally AdSs black brane metric, 

ds 2 = — 2 u^x^dx^dr — r 2 f (b (x a )r) u /1 (x a )u l/ (x a )dx fl dx 1 ' + r 2 A /ll/ (x a )dx fl dx l/ , (3.10) 

where /(?’) = 1 — 1/r 4 . Hence r = 1 is the horizon. The conformal boundary is at r = 00 . 
The velocity held Un(x a ) is normalized in a standard way 

g^(x a ) Utl (x a )u u (x a ) = -1. (3.11) 

The temperature held T(x a ) is related to the parameter b(x“) via 

T(x a ) = * , , (3.12) 

v ' 7rb(x a ) v ' 

which is identihed as temperature of the boundary CFT. If ^(x"), b(x“) and g^, v {x a ) are 
taken constant 

Ufj,(x a ) ---> u^, b(x“) ---> b, gp,v(x a ) ---> g /w , (3.13) 

the line element (3.10) does form a class of solutions to the bulk Einstein equations (3.2). 
As r —> 00 , the metric (3.10) approaches asymptotics 

ds 2 ■ r ~ 5 ~° c ’> — 2 u fl (x a )dx ll dr + r 2 g^ jU (x a )dx >l dx 1 ' + O ^“ 2 ^) > (3-14) 

2 Due to differences in signature convention of the metric tensor, the Riemann curvature of [43] has an 

additional overall minus sign. 
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which identifies gfi V {x a ) as the background metric for the boundary theory. 

To study fluid dynamics we are to consider velocity and temperature as given, but 
at this stage arbitrary functions of boundary coordinates. Then the element (3.10) ceases 
to solve the Einstein equations. In order to find a solution, we have to add to (3.10) a 
correction which is to be determined by solving the Einstein equations (3.2). Our main 
goal is to compute the stress-energy tensor with all the derivative terms linear in the fluid 
dynamical variables and boundary metric perturbation resummed. These variables can be 
expanded to linear order 

u^(x a ) = ^-1 + ^ eh 00 (x a ), etq(x Q )^ + O (e 2 ) , 

g^{x a ) = rj^ + eh lx „(x a ) + G(e 2 ), ( 3 ' 15 ^ 

b(x") = b 0 + ebi(x Q ) + 0(e 2 ), 


where e is an order counting parameter and will be set to one in the end. As seen from (3.15), 
gravitational perturbations can induce fluid flow 3 . The constant bo corresponds to the 
equilibrium temperature of the boundary theory and will be set to one from now on. The 
seed metric, i.e., a linearized version of (3.10) is 


ds z eed = 2 drdv — r z f(r)dv 2 + r 2 dx 2 


— € 


2 Ui(x a )drdx l H—~ Ui(x a )dvdx l -|—ybi (x a )dv 2 


1 


(3.16) 


+ 0(e 2 ). 


+hoo{x a )drdv H— ^ho d (x a )dv z — r z h^(x a )dx fJ 'dx v 

Denoting the metric correction as ds 2 OII , the full metric is then formally written as 

ds 2 = G MN dx M dx N = ds 2 eed + ds 2 orr . (3-17) 


For the metric correction we choose the “background held” gauge [27] 


G rr — 0, G r n oc Tr 




= 0 , 


(3.18) 


where G corresponds to the first line of (3.16) and stands for the metric correction 
ds 2 OII . Then the most general form of ds 2 orr can be cast into the following form 

ds 2 co rr = e ^—3 hdrdv + ^dv 2 + r 2 hdx 2 + jidvdx 1 + r 2 ctijdx 1 dx^ , (3.19) 

where aij is a traceless symmetric tensor of rank two. We would like to stress that all 
the components {h, k , ji, a t j } are explicit functions of the bulk coordinates {x“,r} and, 
through the Einstein equations, also become functionals of Ui(x a ) and h fiu (x a ). 

In order to solve for the metric correction, we have to impose appropriate boundary 
conditions. The relevant ones have been discussed in details in [1, 2], For completeness, 
we briefly summarize them here. The first type is a regularity requirement for all the 

3 Recently, an opposite effect, generation of gravitational waves from sounds was considered in [46, 47]. 
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components {h, k, ji, a^} over the whole range of r. Second, since the boundary theory is 
supposed to live in a fixed spacetime with given metric g liv (x a ). near the boundary r = oo 
we require 

h < O(r 0 ), k < 0(r 4 ), ji < 0(r 4 ), a tj < O(r 0 ). (3.20) 

Finally, the remaining ambiguity is fixed by defining the fluid velocity u fl (x a ) in Landau 
frame, 

u ,l T lw = -eu v =■=> = 0 (3.21) 

which, within the linear approximation, is equivalent to specifying II^jA as transverse 

n ( 0 o> = = 0 ; n (lJ > ^ o. (3.22) 

In Appendix C, the tensor T /iU (r) is explicitly expressed in terms of the metric compo¬ 
nents (3.19). So, once the metric corrections are found from the Einstein equations, it is 
straightforward to derive the stress-energy tensor from (C.1,C.2). 


4 From gravity to fluid dynamics 


In this section, we derive the boundary fluid dynamics through solving the gravitational 
equations in the bulk. The expressions (3.2) are a set of fifteen equations, with one redun¬ 
dant component. The remaining equations can be split into ten dynamical equations and 
four constraints. The dynamical equations are to be solved for h, k, ji and Q!,y, which are 
precisely ten unknown functions in the metric corrections (3.19). Thus obtained solutions 
would still functionally depend on the, yet unspecified, velocity field and background met¬ 
ric. For the boundary fluid, solving for dynamical components of (3.2) would mean that 
we first construct an “off-shell” stress-energy tensor. It is, however, worth emphasizing that 
this procedure is sufficient to uniquely fix the transport coefficient functions. 

The remaining four constraint equations act as constraints in the space of velocity and 
temperature fields, for which solutions to the gravitational equations exist. These equations 
are found to coincide with the Navier-Stokes equations for the resummed hydrodynamics. 
In other words, satisfying the constraints is equivalent to putting the stress-energy tensor 
“on-shell”. We demonstrate this equivalence in Appendix E, where the conservation laws of 
T^ v are shown be to consistent with the constraints. 

4.1 Dynamical Einstein equations 

The first dynamical equation is E rr = 0, which yields 

5d r h + rd 2 h = 0. (4.1) 


The 

The 


asymptotic requirement (3.20) and the fluid frame convention II(oo) = 0 lead to h = 0. 
function k can be found from E rv = 0, 


d r k =2r 2 du - 2r 2 d k h 0 k ~ ^rd 2 h 00 - ^r (didjhij - d 2 h kk ) 

1 2 1 1 

+ r 2 d v h kk + -rd v du - T^dj ~ -rdidjCXij - 7 ^d r dj, 


(4.2) 
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which is coupled to j % and a.ij. Our strategy is the same as that in [2]: we will first solve 
for the functions j, and ay. and then integrate the equation for k. 

The equation for j t is derived from E r i = 0, 

0 = rd'rJt - 3 d r ji + 3r 2 d v iii + r (< d 2 ut - didu ) - ^ r 2 dih 00 + r 3 d r d 3 a i3 , (4.3) 

which is coupled with ay. The equation for the latter is determined from Eij = 0, after 
some massaging, which is presented for the flat space in [2], Extension to a curved space is 
straightforward and here we quote the final equation only 


0 =(r 7 — r 3 )d^aij + (5 r 6 — r 2 )d r a.ij + 2 r 5 d r d v aij + 3 r 4 <9„cty 
+ r 3 (d 2 aij - did k a lk - d 3 d k a lk + 

+ (1 - rd r ) (dijj + djji - ‘^Sijdj'j + ( r 3 d v + 3r 4 ) (diUj + djm - ^ 5ijdu ^ 
+ 3r A d v (hij - ^Sijh kk ^ - r 3 (didjhoo - ^8 l3 d 2 d 2 h 00 
3 r (dihoj T djhoi —8ijd k ho k 

T v ^d~hij did k hjk djD k hj k T 8j0jh kk — 8j 3 0 h kk T ~4y 0 k Dih k j ^ . 


(4.4) 


Eqs. (4.2, 4.3, 4.4) are partial differential equations. Furthermore, they all have source 
terms, which depend on the velocity held m and background metric h^ u . Those are arbitrary 
at this moment. Thus, solutions of (4.2, 4.3, 4.4) are functionals of ry and h^ v . Our strategy 
for solving these equations goes in two steps. First, due to the linear approximation, we 
Fourier transform these equations with respect to the boundary coordinates. In this way 
we cast these equations into ordinary differential equations with respect to the holographic 
coordinate r. As a second step, we rid off the functional dependence on ty and h /w . This is 
done by decomposing the vector j t and tensor aij using base tensor structures constructed 
from Ui and h^ v . These structures are classified according to SO( 3) symmetry and are 
listed in Table 1. We end up with a large system of (partially coupled) ordinary differential 
equations for the coefficient functions. 

To proceed, ji and ay are uniquely decomposed as 


ji — I i u?; + V28idu T V3/10 i T Vidid^hok T Vjjclj/ioo T V§8ih kk T FVSfchj/j T Fs 8id k dih k i, 

a ij = dy + T^^ij + 2T^di 3 + T^ij + T§<pij + T$Xij + 2T7^y + Tgfaj + Tgry-, (4.5) 

where V) and T) are short notations for Vi(oj,qi,r) and Ti(u,qi,r). Since the temperature 
variation tq never appears in any of the dynamical equations, the tensor structures con¬ 
structed from bi in Table 1 do not enter the decomposition (4.5). Substituting (4.5) into 
eqs. (4.3, 4. 4), we arrive at the ordinary differential equations for V) and Tj. In what follows, 
we present the equations grouped into decoupled sectors. 
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Scalar 

Vector 

Tensor 

bi 

Ui 

Tij = hij g ^ij hkk 

hoo 

hoi 

&ij = \ (diUj + djm - | dijdu) 

hkk 

<9jbi 

hij = 2 ( pihftj ~\~ djhoi 3$ijdkhok) 

du 

dihoo 

TUij = didjhi - \6ijd 2 hi 

dkhok 

di^kk 

Lpij = did j hoo 3 $ijd hoo 

Oj hij 

didu 

Xij — did j hkk 3 $ij d hkk 


didkhok 

^Pij — 2 {didkhjk H - djdkhik 3 $ijdkdihki ) 


didkdihki 

7tij = didjdu — ^Sijd 2 du 

TTij = did j dkhok 3 $ijd dkhok 

(fiij = didjdkdihki 3 &ijd dkd\hki 


Table 1. Up to linear order 0(e), we exhaustively list all the basic tensor structures constructed 
from the fluid dynamic variables bi and classified according to the S'0(3) symmetry. 

I: {V u V 2 , T,, T 2 } 

' 0 =rd 2 Vi - 3d r Ui - 3iur 2 - q 2 r - q 2 r 3 d r T l , 

0 =rd 2 V 2 - 3d r V 2 - r + ^ r 3 d r Ti - | q 2 r 3 d r T 2 , 

o o 

0 =(r 7 - r 3 )d 2 T\ + (5r 6 - r 2 )d r Ti - 2iur b d r Ti 

— 3iu:r' i Ti + Vi — rd r V\ — zwr 3 + 3?’ 4 , 

0 =(r 7 - r 3 )d 2 T 2 + (5r 6 - r 2 )d r T 2 - 2 iur 5 d r T 2 

- 3iur 4 T 2 + 2U 2 - 2rd r V 2 + \q 2 r 3 T 2 - ^r 3 Ti. 

3 3 

This first sector decouples from the metric perturbation and was explored in [1, 2]. The 
holographic RG-flow type equations (4.6) determine the viscosity functions r] and 
The remaining sectors are new. They all emerge due to the metric perturbation. 

II: {V 3 , V 4 , T 3 , T 4 j 

' 0 =rd 2 V 3 - 38 r V 3 - q 2 r 3 d r T 3 , 

0 =rd 2 V 4 - 3d r V 4 + * r 3 d r T 3 - ^ q 2 r 3 d r T 4 , 

0 =(r 7 - r 3 )d 2 T 3 + (5r 6 - r 2 )d r T 3 - 2iur 5 d r T 3 

- 3iur 4 T 3 + U 3 - rd r V 3 - 3r 4 , ' ‘ ’ J 

0 =(r 7 - r 3 )d 2 T 4 + (5r 6 - r 2 )d r T 4 - 2iur 5 d r T 4 

— 3iur 4 T 4 + 2V 4 — 2rd r V 4 + -g 2 r 3 T 4 — -r 3 T 3 . 

3 3 
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HI: {V 5 , T 5 } 


0 =rd 2 r V\ 5 - 3<9 r V 5 - j|r 2 - q 2 r 3 d r T 5 , 

< 0 =(r 7 - r 3 )S 2 T 5 + (5r 6 - r 2 )d r T 5 - 2iur 5 d r T 5 (4.8) 

- 3 iojr l T$ + \q 2 r 3 T 5 + 2V 5 - 2rd r V 5 - r 3 . 

O 

IV: {T 9 } 

0 =(r 7 - r 3 )d 2 T 9 + (5r 6 - r 2 )O r T 9 - 2iujr 5 d r T 9 

— 3iur 4 T 9 — q 2 r 3 T 9 — 3 zwr 4 — q 2 r 3 . ^ ^ 

V: {V 6 ,T 6 } 

0 =rd 2 V 6 - 3d r V 6 - ^r 3 d r T 9 - ^ q 2 r 3 d r T 6 , 

< 0 =(r 7 - r 3 )d 2 T 6 + (5r 6 - r 2 )d r T 6 - 2iojr 5 d r T 6 (4.10) 

- 3zwr 4 T 6 + 2V 6 - 2rd r V 6 + - q 2 r 3 T 6 + ^r 3 T 9 + r 3 . 

00 

VI: {V 7 , T 7 } 

' 0 =rd 2 V 7 - 3d r V 7 + r 3 d r T 9 - q 2 r 3 d r T 7 , 

< 0 =(r 7 - r 3 )d r T 7 + (5r 6 - r 2 )d r T 7 - 2iivr b d r T 7 (4.11) 

- 3iur 4 T 7 + V 7 - rd r V 7 - r 3 T 9 - r 3 . 

VII: {V 8 , Tg} 

0 =rd 2 V 8 - 3d r V 8 + ^r 3 d r T 7 - q 2 r 3 d r T s , 

< 0 =(r 7 - r 3 )d 2 T 8 + (5r 6 - r 2 )d r T 8 - 2iur 5 d r T 8 (4.12) 

- 3iur 4 T 8 + 2Vg - 2rd r V 8 + (m 2 r 3 r 8 - ^r 3 T 7 . 

OO 

The transport coefficients are completely fixed by integrating these equations from the hori¬ 
zon to the boundary. We, however, postpone this integration step until subsection 4.3. In 
the next subsection 4.2, we focus on the fluid stress-energy tensor, which so far is expressed 
in terms of h, k, ji and ctij (see (C.1,C.2)). Via the decomposition (4.5), the stress tensor 
can be written in terms of momenta- and r-dependent coefficient functions V) and T) mul¬ 
tiplying the tensorial structures listed in Table 1. Furthermore, we are mostly interested 
in these coefficient functions near the boundary. Eventually we will be able to map certain 
pre-asymptotic behavior of these coefficient functions onto momenta-dependent transport 
coefficients. To this goal, we first analyze the above equations near r = 00 . 

4.2 Stress-energy tensor from near the boundary analysis 

We are to study the large r behavior of the metric corrections, which makes it possible 
to cast the stress tensor into the form ( 1 . 1 , 1 . 8 ), but with the transport coefficients left 
undetermined. The latter cannot be determined from asymptotic considerations only, but 
require a full integration over the holographic coordinate. 
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We begin with the first sector {Vi, T\, V 2 , T 2 }, which was already analyzed in [2]: 


V 1 -iur 3 + 0[- 

r 


rp'-t \ rpO 


V 2 


r—>■ 00, 1 2 1 /r\ ( 1 \ m r—>■ oo t ^2 ^ / 1 


t - 3 r+ °Vh 


The momenta-dependent coefficients fi and £2 are related to the viscosities [1, 2] 

V = —4ti, C = -4*2- 

In the limit r —> oo, the functions {V3, V4, T3, X4} have the following expansion 

V 3 —> + \iuq 2 r + T779 2 (g 2 - w 2 ) logr + -*-4(w 2 + 3<? 2 ) + C 

4 o Id 54 

1 iu q 2 — u 2 t-t 1 . 9 9 . , ^ /logr 

^ H - I6H '° g r + ° ( 

1 


T 3 

V 4 

t 4 


\ r 
logr 


- 9 1 l.o o. a; 2 + 3q 2 , 

l 2 r -I8'“ r -l92 ( ‘" -" ) + ^— l0gr + °, r 

1 t 4 iuj logr „„ /logr 

+ 4+ _ , +0 1 


6 r 3 r 


12r 4 


As r -> 00, all the other functions have the asymptotic behaviors of the form, 

T r 1 3 . 1 2 1 ■ 2 1 1 • 21 |/-o ( l°g r 

V5 —t - -r + -q r - — iuq + —logr + O I 

T 5 ^4^ + l 2 ^ + ^ + 48^ ((? " 3a; )l0gr + °V^ _ 

I/ 6 —> - + 4(^2 _ w 2 )r _ _L^( g 2 - u 2 ) 

+ ^M <? 2 - ^ 2 ) log r + O 

T t 1 ■ 2 , 1 2 1 ■ , 2 1 o 2\ 1 • , 2 2 \ 1 ,,-0 /logr 

Vi —>-nor + -lo r — —iuj(uj + 3 q ) — — iiu(q — u Jlogr + O 


6 


64 


16 


^8 

Tq 

t 7 

t 8 

t 9 


1 5 . 1 /logr 

► — r- iu --zculogr + Cl - 

18 96 24 B \ r 

1 IUJ 1 tfi 1 , 2 2 m , m 4°S r 
’4r 2 12r 3 + r 4 48r 4 ^ U ^ logr + ° \ r 5 

1 . . *7 1 , 2 2 m . M|4°g r 

^"^ + ^3 + ^4 -^)logr + G 

t 8 log r- / log r 

V 4 + 24r 4 + 

‘ “ y “ ^2 (<? +w) + ^3*w(3g - W ) + ^4 

1 , 2 2\2 1 m 4 °S r 

+ I 6 ^ (C? } logr + 0 lT^ 


(4.13) 


(4.14) 


(4.15) 


(4.16) 
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Then, from equation (4.2), we can obtain the large r limit of k 

k r -°°> | (r 3 + iur 2 ) du - ^r 3 d k h 0k + * r 2 d 2 h 00 - ^r 2 (didjhij - d 2 h kk ) 

+ xv [-4 iu>q 2 d k h 0k - 2 q 2 d 2 h 00 + (u 2 - q 2 )d 2 h kk 

5 (4-17) 

+(q 2 - 3w 2 )9fc^h fc /] logr - — [4*wg 2 3 fc /i 0 fc + 2q 2 d 2 h 00 + 

(q 2 - u 2 ) d 2 h kk + (3a ; 2 - q 2 )d k dih k i] + O . 

The integration constants in V)’s and k are fixed by the convention Il^gx = 0 and n^o) = 
0, respectively. The asymptotic requirements (3.20) were used to determine the leading 
behaviors in the near-boundary expansion. Our problem is now reduced to finding a set 
of pre-asymptotic coefficient functions ij, which propagate into final expression for the 
boundary stress tensor (eq. (4.18)). These coefficients can be computed only through full 
integration from the horizon to the boundary. The regularity conditions at the horizon 
(r = 1 ) provide sufficient initial data to uniquely fix all the ti s, and consequently all the 
transport coefficient functions introduced in ( 1 . 8 ). 

We also notice the logarithmic branches in (4.15,4.16,4.17), which emerge solely due 
to the boundary metric perturbation h^, as is clear from (4.5). The logarithms appear 
starting from the fourth order in the boundary metric derivatives, in full agreement with 
computations of [43, 44, 48]. 

The stress-energy tensor for the boundary CFT is obtained by substituting the above 
large r behaviors into (C.1,C.2), 


Too = 3 — 12ebi — 3e/ioo, 

T 0i = T i0 = ~4eui + eh 0i , 

Tij = 5ij( 1 - 4ebi) + ehij + e {8ti<nj + 4f 2 7%■ 


+ 

+ 

+ 

+ 


8(ii + t 3 ) + —iuj(3q 2 - 7a; 2 ) 


dij + 


4 (£2 + £ 4 ) + — 
36 


4i5 “ 144 ^ 2 + 21uj2) 


<fij + 


7%' 
2 , , 2 ’ 


8t? — ^(w 2 + 3 q 2 ) 


ipij + 4f 8 - — 


2 iut 2 + 4t 6 + — (13<T ~ u ) 


7 


Xij 


1 


72 ' T%3 
2 2 


Aiuoti + Atg + -(-(7a; 4 — 9 q 4 — 6 u 2 q 2 ) 
48 


T, 




The tensors cr,j and 7 are covariant versions of those defined in [1, 2], 


(Tij — I + V jUi gij\7u j — (Jjj dij iiOTij, 


iiij = VjVjVrt - -gijV 2 Vu = itij - - -iuxij, 


(4.18) 


(4.19) 


where the linearization approximation (3.15) was employed in expressing them in terms of 
the basic structures defined in Table 1. Combining non-derivative parts in (4.18), we arrive 
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at the standard expression for conformal ideal fluid, 


^ eal = ^(^ + 4 V ,). (4.20) 

The expansion coefficients ti s are not fully independent. There are four constraints among 
them, derivation of which is deposited to Appendix D: 

0 = tg — iujt 3 — q 2 t 7 , 

0 = 2 (ii + t 3 ) - 2 iut 5 - q 2 {t 2 +U), 

„ (4.21) 

0 = 2 (t 6 + tj) — iujt^ — 2 q 2 t%, 

0 = (*5 - h) - iuU - q 2 t 8 . 

Eq. (4.18) can be rewritten in terms of the Weyl tensor C^ av p and its derivatives. With the 
help of Appendix A and after some algebraic manipulations including the constraints (4.21), 
the dissipation tensor takes the form ( 1 . 8 ) with the transport coefficients r] etc being 
expressed as combinations of the coefficients t t 's. 

rj = -4ti, C = —4*2, 0 = — 6^2 

K = \ [—48(^5 + f 6 ) + 24icu(t 2 + 2 t 4 ) + w 2 (48 1$ - 1 ) - q 2 ] , 

o , 

4 (4-22) 

p = - - [12 (t 2 + t 4 ) + iu{ 1 - 24t 8 )] , 

Z = ^(7-288t s ). 

Notice that, in addition to the explicit u> and q dependencies present in (4.22), all the 
coefficient functions ti also depend on these momenta. 

As has been emphasized above, the stress-energy tensor (4.18) is off-shell: so far we have 
managed to establish a map between the dynamical components of (3.2) and the transport 
coefficients. In our procedure, we have specified neither the velocity nor the boundary 
metric perturbation. We expect that the conservation laws \/^T^ u = 0 should emerge from 
the constraints in (3.2), more precisely, from the constraints E vv = E v i = 0. We have 
checked that it is indeed the case. The relevant computations are deferred to Appendix E. 

4.3 Gravitational susceptibilities of the fluid 

We have now reached the point where we are to solve the bulk equations (4.6-4.12). These 
differential equations have to be evolved from the horizon to the conformal boundary with 
appropriate boundary conditions imposed at both ends. For the sector (4.6) this calcu¬ 
lation has been performed in [2] and revealed the viscosity functions. Here we apply the 
methods developed in [2] to other sectors. We first perform analytical analysis in the hy¬ 
drodynamic limit and then full numerical study for generic values of momenta. Below we 
merely summarize the results, while referring the reader to [ 2 ] for calculational details. 
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4.3.1 Analytical results: hydrodynamic expansion 

We start with the hydrodynamic limit ui, qi <C 1 and solve the equations (4.6-4.12) pertur- 
batively. We introduce a formal parameter A by u —»• Xu, qi —> Xqi, and expand Vi’s and 
Ti s in powers of A 

OO OO 

Vi{u,qi,r) = ^r\ n v} n \u,qi,r), Ti(u,qi,r) = ^A n T- n) (u, q t ,r), (4.23) 

A =0 A =0 

where A will be eventually set to unity. At every order in A, a system of ordinary second 
order differential equations can be integrated [2], The expansion coefficients V- n ' > and T^ n ' > 
are expressible as double integrals, which are summarized in Appendix F. Their r —> oo 
asymptotic behavior reveals the coefficients tj’s: 

t i = ~\ + —^+ -^2 {6 q 2 + u 2 [6vr- 7T 2 + 12(2 - 3 log 2 + log 2 2)] } H-, 

1 , , 1 log 2 - 1 

t 2 =-^ (5 - 7T - 21og2) H-, t 3 = - --- iu^ -, 

U = (1 + vrT21og2) H-, t 5 = ^ (4 + vr — 2 log 2) iu H-, (4.24) 

*6 = ^ (10 + 7T - 10 log 2) iw H-, t 7 = ^ ^ (4 log 2 - 3)iu H-, 

~ -288- 1 -’ td ~ 4 lUJ + g + ~ l) w ^] -• 

When substituted in (4.22), this leads to the expansions (1.9) for the viscosities and GSFs. 


4.3.2 Numerical results: all-order resummed hydrodynamics 


We are now to solve the equations (4.6-4.12) for generic values of u and q 2 . Since the 
boundary conditions are imposed at different points, we resort to the shooting technique 
and solve them numerically. As in [2], we first start with a trial solution (initial condition) 
at one boundary (horizon) and integrate these equations to the second boundary (the 
conformal boundary). Thus generated solution should match the boundary conditions at 
the end of the integration. So, the trial solution has to be finely tuned in order to have the 
boundary conditions at the end of the integration satisfied. This fine-tuning procedure is 
mapped into an optimization problem. 

The viscosities r/ and £, originally computed in [2], are shown in figures 1 and 2. r/ and 
£ are observed to vanish at very large momenta. As has been pointed out in section 2 , it 
is worth looking at the corresponding memory function defined through the inverse Fourier 
transformation, 


fj ( t,q 2 ) 



(4.25) 


The viscosity r) is weakly dependent on spatial momenta, meaning it can be thought of as 
almost a local function of space coordinates. For this reason, we have chosen not to Fourier 
transform in q and present the memory function in a mixed ( t , q 2 ) representation. In 
figure 3, we show the time dependence of the memory function fj(t,q 2 ) and compare it with 
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Figure 3. Memory function rj ( t,q 2 ). Left: 3D plot as a function of time t and q 2 . Right: 2D plots 
as functions of time f; the dashed curve corresponds to Israel-Stewart hydrodynamics with 770 = 1 
and the relaxation time tr = (2 — log 2) /2; the solid curves display our result at different q 2 (from 
the rightmost, q 2 = 0,1, 2, 3). 


the Israel-Stewart hydrodynamics, for which the viscosity function 77/5(0;) = 770/(1 — iwi/j) 


19 - 






Figure 5. GSF p as function of ui and q 2 . 



Figure 6. GSF £ as function of u> and q 2 . 


and the memory function is 


Visit) 


^%oe- t/T *0(f). 


tr 


(4.26) 


Notice that the memory function r) vanishes exactly for negative times, in accordance with 
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the causality requirement, as discussed in section 2. For positive times, the memory function 
decreases with time, followed by damped oscillations, and eventually approaches zero at very 
late times. 

We are now to present the results for the GSFs. As has been proven above, 9 = 3£/2. 
Numerical results for the remaining GSFs are displayed in figures 4, 5 and 6. Compared 
to the viscosities, we observe a noticeably different behavior of the GSFs as functions of 
momenta. Particularly, some components do not vanish at large momenta. Below we shall 
provide interpretation of such behavior. 

k is found to depend weakly on spatial momentum qi, somewhat similarly to the viscosi¬ 
ties 7] and £. Higher order derivative terms reduce Re [k] until the point (uj ~ 2.6, q 2 = 0}, 
where Re[«;] reaches its minimum. At larger frequencies, Re [k] oscillates around zero and 
then it starts to grow around uj ~ 5. The imaginary part of k, increases from zero to its 
maximum at {cu ~ 1.8, q 2 = 0}. Then, Im[fc] decays very fast until {uj ~ 4.5, 9 = 0} where 
it starts to oscillate around zero. At very large momenta, the real part of k keeps growing 
while its imaginary part vanishes. 

The GSFs p and £ behave rather differently from their cousin k. When uj < 2, Re[p] 
drops sharply as the spatial momentum is increased. For larger uj, the real part of p first 
decays rather slowly until the point {w ~ 2.0, q 2 = 0}, where Re [p] starts to drop quickly 
and eventually approaches zero at very large momenta. In contrast, the imaginary part 
of p grows fast at large frequencies. The spatial derivatives also affect Re[£], in the same 
region of uj < 2. With uj increased, Re[£] grows fast till uj ~ 4, where it tends to approach 
a constant value of 1.2. The imaginary part of £ weakly depends on spatial momentum. It 
first decreases from zero till certain minimum around {w ~ 3.2,9 = 0}, where it starts to 
increase as function of uj but approaches zero starting from uj ~ 5. 

We do not have a good intuitive explanation of the observed behaviors of the GSFs 
from small to mid momenta uj ~ q ~ 5. Nevertheless, the large momentum asymptotic 
region can be understood. In the limit uj, q —> oo we expect to recover zero temperature 
limit of the theory. The viscosity functions vanish in this limit, whereas the GSFs do not. 
While the role of each individual GSF coefficient is not obvious, their combinations enter 
as residues in the expression for two point correlators [3]. We review the correlators in 
Appendix E. Thus we expect the GSFs to play a role in recovering the zero temperature 
limit of the correlators, which are known analytically and also quoted in Appendix E. We 
have checked numerically that, starting from uj ~ q ~ 5, the correlators computed using 
the GSFs (eqs. (E.6,E.9,E.12)) reproduce well their zero temperature counterparts (E.13). 

From the zero temperature limit of the correlators, we can analytically deduce the 
asymptotic behaviors of the GSFs, 


k ~ (a; 2 - 9 2 ) 2 log ( uj 2 - q 2 ), 

p^- (uj 2 - 9 2 ) log (uj 2 - 9 2 ), (4.27) 

UJ 

£, ~log (oj 2 - 9 2 ) 

which agrees with the numerical results presented above. 
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5 Summary and discussion 


In this work we introduced (non-dynamical) metric perturbations of the flat boundary, and 
following the the setup of [1, 2] obtained an asymptotically locally AdS ,5 black brane solution 
to the Einstein gravity. Our construction is limited to linear order in amplitudes of the fluid 
velocity and boundary metric perturbation, but exact in terms of the derivative expansion. 
This gravity solution is holographically dual to a relativistic conformal hydrodynamics in 
a weakly curved background spacetime. This dual hydrodynamics is equivalent to full 
linear response theory with the background metric perturbation playing the role of external 
disturbance. Since all order derivatives are included in the stress tensor, the effective theory 
constructed here is causal and expected to be free of any instabilities. 

The constitutive relation (1.8) for the stress-energy tensor emerged from the dynami¬ 
cal components of the bulk Einstein equations, with holographically determined values of 
transport coefficients. Overall we introduced six transport coefficient functions: two vis¬ 
cosities and four GSFs, transport coefficient functions associated with ffuid’s response to 
curvature perturbations. These terms vanish in fiat space and do not affect hydrodynamic 
equations. When a non-zero gravitational perturbation h^ u is turned on, one can consider 
it as an external force such that the dual fluid becomes forced in Minkowski space [30]. The 
arguments of [30] indicate that different gravitational forces are able to stir the fluid into 
flows with velocity differences of order one. 

It appears that the GSFs have a larger role beyond being naive linear response coeffi¬ 
cients. As was first noticed in [40] and further elaborated in [3], GSFs enter the expressions 
for thermal correlators of T^ u . In Appendix E, the correlators are derived using the linear 
response theory from the constitutive relation (1.8). The shear and sound wave modes pre¬ 
dicted from hydrodynamic equations in flat spacetime appear as poles in these correlators. 
Residues of the thermal correlators get contributions from both the viscosity r] and all the 
GSFs (see (E.6,E.9,E.12)). We have found that some of the GSFs do not vanish at large 
momenta. This is consistent with the expectations put forward in [3]. It was argued there 
that zero temperature expressions for the correlators have to be recovered from the thermal 
ones in the large momenta limit and that is what the GSFs take care of. The thermal 
correlators incorporate information about both thermal physics related to matter flow and 
the vacuum contribution due to pair production. Furthermore, there are interference effects 
between these two types of physics, as encoded in various GSF terms. 

A remark about Weyl symmetry is in order. The stress tensor of conformal fluid 
is supposed to transform covariantly under the Weyl rescaling —> e -2 ^^ [29, 30]. 

In [40], the symmetry was used to impose selection rules on the form of second order 
relativistic conformal fluid. It was subsequently exploited in [3] in order to constrain the 
general form of linearly resummed constitutive relation 1 . Both the Navier-Stokes term 
and the Weyl tensor transform as needed under the Weyl transformation. 

It is, however, obvious that applying extra derivatives on them breaks the symmetry. It 
turns out that the linearization does not commutate with the Weyl transformation: the 
nonlinear terms like A\Qy... (see Appendix B for introduction to Weyl-covariant formalism) 

4 Thanks to linearization, the Weyl anomaly is irrelevant for the discussion. 
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do generate some linear pieces in the Weyl transformed spacetime g /1L/ = e 2< ^g^ u - The Weyl- 


covariance of the fluid stress-energy tensor is to be recovered via a non-linear completion 


of the expression (1.8) [3]. In order to have a manifestly Weyl-covariant form of the stress- 


tensor, we can substitute the normal derivative by Weyl-covariant long-derivative 
introduced in [49] (see Appendix B). With this procedure we would also recover a significant 
fraction of the non-linear terms neglected by our analysis. 


We have been emphasizing the relation between all-order gradient resummation and 


causality. Particularly, in section 2 we argued that causality induces constraints on posi¬ 
tions of the poles of the transport coefficients. We demonstrated this effect by explicitly 
computing the viscosity memory function, which indeed does not have any support in fu¬ 
ture. In a forthcoming publication [53], we study all order linearized fluid dynamics whose 
dual gravity involves a Gauss-Bonnet correction. We find that the value of the Gauss- 
Bonnet coupling can be constrained by causality requirement on the linearly resummed 
fluid dynamics. 

We hope our results may share some generic features with those of realistic QCD and 
could be helpful in understanding the nature of hydrodynamic gradient expansion. In 
fact, when resumming all order derivative terms, we explicitly go beyond the conventional 
hydrodynamic limit and account for non-hydrodynamic modes in our construction. At 
momenta of order one the hydrodynamic and non-hydrodynamic modes mix and only full 
resummation makes sense. 

The viscosity memory function computed here could be useful for phenomenological 
applications. We have illustrated how different it is from the naive second order ansatz a la 
Israel-Stewart. In [50], higher order gradients were argued to have large phenomenological 
effects. So, it is quite important to better explore these effects by simulating hydro evolution 
of quark-gluon plasma with our viscosity memory function, despite its obviously limited 
applicability to real QCD. 

A Boundary curvatures 

In this appendix, we summarize expressions for boundary curvatures that are useful in 
constructing T^ u . We work with the mostly plus signature for the metric tensors. Thus, 
the Riemann curvature of g^ u is conventionally defined as 



(A.l) 


where T^, is the Christoffel symbol. To linear order in R^ v and R have the expressions, 



Rq i — 2 ( ^ T didkhok T d v dkhik dv^ihkk) 
— 2 [ <9 hij T (didkhjk T dj^khik) T didjh oo 



(A.2) 


d v (dihoj T djhoi) + d v hij didjhkk ] , 

R d hoo & ^kk T &k&ihki ‘^dydkh'Ok T & v hkk 
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which were used to compute the logarithmically divergent piece The Weyl tensor is 
defined as 

C\ fieri/ — R\furu 2 (^A crdfi.i' R\v9fJ.cr Rfiad\v RfiudXa) T {.9\cr9fiv 9\i'9[J.cr') ■ (A.3) 

In terms of the basic structures listed in Table 1, the tensors introduced in (1.8) are 

Uj) = CTij (7jj + V(jV U = 7Tjj 7Tjj T T^pvXiji 

u a uPC( iaj )i 3 = - * (3 2 Tij + - 2^ + + Xij + 2 d v aij) , 


u a V /3 C( ia j' )/3 — - [-2d v Ciojo + (C i0 jk + Cikjo)] 

— ~ (d v Tij d v d Tij T dy'ipij + d v (pij 2d v dij T d Cij 7r^) , 

V a \ / ^C^ ia j' ) i 3 = — 9,;3 fc (Ciojfc + Cj fc: ,o) + dkdiC( ik j)i 

= ^ (6 dld 2 Tij - - 33 4 rij - 69 2 ^y - 3<9 2 (^ + + Gc^d,: 


(A.4) 


-6 d v d 2 aij + 4:d v Ttij + d 2 ^ - d 2 Xij + 69 2 ^i - 2^) , 


^ ’o.Tt{ij) — „ ( d v T~ij d v d Tij + 2d v ipij + d v if>ij 2d v dij d v Xij) j 


where the linearization approximation (3.15) has been applied. 


B Weyl-covariant formalism for relativistic conformal fluid 


In this appendix, we briefly summarize the Weyl-covariant formalism of [49]. The main 
ingredient is the Weyl-covariant derivative operator T>\. Consider a Weyl-covariant tensor 
Qy... which transforms homogenously under the Weyl rescaling, 

Qu..'. = e~ w<t> Q%;:. under g, iu = e 2(t> g^ v . (B.l) 


The Weyl-covariant derivative is defined as 


'DxQ^:: =v A q$:\ + wA x Qt::: + [< gxa A? - s»a q 

- [g Xu A a - 5 x Av - 8$A x ] Q£: 
where A^ is the Weyl-covariant connection and is dehned as 

Afj = U V a^fi 3^V„« ■ 


Qu.:: + • • • 


Under the Weyl rescaling, A/, transforms in the same way as 17(1) connection, 


(B.2) 


(B.3) 


Aft = A,, + d^, under = e 2<l> g llu . (B.4) 

Notice that under linearization, A ~ d^u an< 4 hence any term which is being multiplied by 
A is formally non-linear. The Weyl-covariant derivative T> x is metric compatible V x g^ = 0. 
The most important feature of V x is that it commutes with the Weyl transformation, 


V x Qi\:: = e- w *V x Qt:: under g^ = e 2 % v . (B.5) 
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C T /iU (r) in terms of metric corrections 

In this appendix, we provide the expression of T /JU (r). In terms of the metric corrections, 
the first piece T£ v (t) is 


Too = 3 “ 12eb i “ 3e/l oo + e Q r 2 didjhij - ^r 2 d 2 h kk - r 2 d 2 h + ^didjOtij 
—9 r 4 h + 3 k — 2 r 3 du + 2r 3 8 k h 0k — r 3 d v h kk + ~dj — 3 r 3 d v h — 3r b d r h S j 


(1 „ 1 

T&j = — 4ertj + e/ioi + e I —r <9j/ioo H— d{k + 4j, 

\2 r 

d h,Qi + —f did k ho k -T —i"* d v d k hi k 
+ i^didj + ^r 2 d v d k a ik ~ r 2 d v dihj , 


3 

- r«9 r ji - r 3 d v Ui - -r 3 dih 
—r d v dihkk ‘2/'~ ^ 


Tjj —dij (1 — 4ebi) + ehij + eSij 

1 


~^r 2 d 2 h 00 + * r 2 d 2 h k k ~ \r 2 d k dih k i 


+r 2 d v d k h 0k - ^r 2 d 2 h kk + \r 2 d 2 h - -^d 2 k - ^—d k dia ki + \d v dj 
2 2 2r A 2 r z 

—r 2 d 2 h + 9?’ 4 /i + k + 2r 3 du — 2r 3 d k ho k + r 3 d v h kk - dj — r 3 d v h 

r 


(C.l) 


+ -d v k + 2 r 5 d r h — rd r k 
r 


+ € 


^r 2 didjhoo - *r 2 d 2 hij - \r 2 didjh kk 


T 2 f (did k hj k T djd k hik) d v (<9j/ioj -\~ djhoi ) + v d v hij 


1 


1 


1 „2o2„, , 1 2 


rdidjh + -^didjk - -r 2 d 2 a^ + -r ( did k a jk + djd k a ik ) 


-^2 dv + d oJi) + 2 r d v a ij - r (d iUj + djUi ) + r 3 (dj/i 0 j + djh 0i ) 


—r 3 d v hij + - (dijj + djji) - r 3 d v ctij - r 3 d, 


r&ij 


where it is important to notice that the terms that explicitly vanish at r = oo have been 
ignored. We also need the expression for T“ u (r) 


na O o2. , q2o2t a 4. oa2a o, , a 2 ■ 


rriCL _ 

J oo —" 


12 


-e {-2d hoo + d v d h kk - d h kk - 3d v d k dih k i + d d k dih k i 
+4d v d 2 d k h 0k ) , 


TSi {3 {did 2 - <9 4 ) h 0i - 2d v d 2 dih 00 + {d 2 + 3d 2 ) did k h ok 

< + {d 3 - d v d 2 ) dih kk - 2d v did k dih k i - 3 {d 3 - d v d 2 ) d k h ik } , ( C ‘ 2 ) 

Tij ( 3 (^ 2 - 3 <9 2 ) Vij + 3 {dv ~ d 2 ) Xij ~ fyij + l2d v Trij 

+185^ (d 2 - d 2 ) dij + 18 {d 2 - d 2 ) - 9 (d 2 - d 2 y T tj 

+5ij [2d 2 h 00 - (d 2 - d 2 ) d 2 h kk + (d 2 - 3d 2 ) d k d,h kl ] } . 

Definition of the tensors etc. can be found in the Table 1. 
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D Derivation of constraints (4.21) 


In this appendix, we derive the constraints (4.21) by suitably combining the differential 
equations of subsection 4.1. Consider the variables 


X\ —Tg — iojT 3 — q 2 Tj, Y\ = —iuV 3 — q 2 V^j , 

X 2 = 2 (Ti + T 3 ) - 2iuT 5 - q 2 (T 2 + T 4 ), Y 2 = 2 {V 4 + V3) - 4ivV s - 2 q 2 (V 2 + 14) , 
X 3 = 2 (Tq + Tj) — iujT. 4 — 2q 2 T$, Y 3 = 214 + 14 — iwl4 — 2q 2 V&, (D.l) 

A4 = T 5 — Tq — iuT ,4 — g 2 T 8 , I 4 = 14 — 14 — iw!4 — g 2 14 + -r 3 . 

These new variables obey similar equations as Tj’s and V^’s, 


0 = (r 7 — r 3 )d 2 Xi + (5 r 6 — r 2 )0 r Xi — 2iur b d r X\ — 3zwr 4 lti + 14 — r<9 r Yi, 

0 = rd^Yx - :id r Y\ - q 2 r 3 d r X 1 . ' ^ ” J 

0 = (r 7 - r 3 )d 2 X 2 + (5r 6 - r 2 )d r X 2 - 2iur 5 d r X 2 - 3iwr 4 X 2 + \q 2 r 3 X 2 

< + y 2 - rd r Y 2 , (D.3) 

0 =rd 2 Y 2 - 3 8 r Y 2 - \q 2 r 3 d r X 2 . 

0 = (r 7 - r 3 )d 2 X 3 + (5r 6 - r 2 )d r X 3 - 2 iur 5 d r X 3 - 3 iur 4 X 3 + \q 2 r 3 X 3 

O 

< + 2 (I3 - rd r Y 3 ), (D.4) 

0 =rd 2 Y 3 - 3 d r Y 3 - \q 2 r 3 d r X 3 . 

O 

0 = (r 7 — r 3 )d 2 X ,1 + (5?’ 6 — r 2 )d r X 4 — 2iur 5 d r X 4 — 3icur 4 X 4 + ^-q 2 r 3 X 4 

O 

< + 2 (Y 4 - rd r Y 4 ), (D.5) 

0 = r9 2 l 4 — 3<9 r I4 — - q 2 r 3 d r X 4 . 

3 


Notice that the equations of JQ’s and Y^s are homogeneous. Under the boundary conditions 
summarized in section 3, the functions IQ’s and Yi s have only trivial solutions just like the 
function h. Furthermore, the identity Yi = 0 does not give non-trivial constraint relation. 
Therefore, we arrive at 

X\ = X 2 = X 3 = X 4 = 0, (D.6) 


which result in the constraint relations (4.21). 

E Navier-Stokes equations and correlation functions 

In this appendix, we first show that the conservation laws = 0 emerge from the 

constraint components of (3.2). We find it more convenient to study certain mixtures of the 
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constraints with the dynamical equations. Particularly, the combination r 2 f(r)E vr +E vv = 
0 is 


11 1 1 r 3 1 

<9 v bi =du — -d k h ok + — 9 2 bi + —— 0 2 /i O o — d 2 hoo + 77 d v h kk 
3 3 Sr 24r 24 6 

- + Y^ r3 ) dvdu + ^dvdkhok - Y^r 3 d 2 h kk - ^dj 

— T^~9 2 k + \d v k + ~!~d v dj — t'- (\ — r^) d r dj , 


12r 


6 r 


12 


(E.l) 


whose large r limit yields 


t^tibi 0 du O d k ho k E „d v h kk . 

3 3 fa 

The combination r 2 f(r)E r i + E v i = 0 leads to 

<9;bi =d v Ui + ^r 3 (didu - d 2 Ui ) + (d 2 h 0i - didjh 0 j) - ^dih 00 - ^r 3 d v dih 00 
E ^r 3 ( d v dih kk - d v djhij ) + ^ r 3 d 2 Ui + ( d 2 ji - 0j0j) + ^d t k - d v j { 

- | r 3 d v djatij E |(r — r 5 )d r djaij - ^ rd r dik + ^r9 r 5„ji. 


(E.2) 


(E.3) 


Taking the large r limit of (E.3), we arrive at 


8ibi - \ 9 ‘ h «> +;h (8i8 “ +382,,i)+ h^ Biau 

E ( 192t s + 3<9„<9 2 - 70 3 ) 5 2 h 0i 

+ ^ (576t 3 + 1152t 4 0 2 - 210 3 - 310, d 2 ) did k h ok 

E ^Y + 219 2 + 0 2 ) d 2 dihoo (E-4) 

+ (192*9 + 192* 7 0 2 + 70 $ - 30 2 0 2 ) 0 fc h ifc 

+ y^Y ( 576t7 + 1152 *80 2 - 0 2 + 90 2 ) did k dih k i 

E yyyy [ 1152i 6<9 2 - 576*9 + ( d 2 - 0 2 ) (210 2 + 0 2 )] dih kk . 

It can be checked that (E.2,E.4) are equivalent to the conservation laws = 0 with T^ v 

given by (4.18). Equations (E.2,E.4) are usually referred to as the Navier-Stokes equations 
for hydrodynamics, which determine the time evolution of the temperature and velocity 
fields once appropriate initial data is specified. 

From the hydro-like constitutive relation (4.18), it is straightforward to construct two- 
point retarded Green’s functions for T^ v , which characterize the linear response of the fluid 
system slightly disturbed from its equilibrium state. In parallel with [3], we present them 
by different channels, as classified according to the 50(3) symmetry. 
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Scalar channel To study the correlation function in the scalar channel, we just need to 
turn on the perturbation h xy (v,z) and particularly the fluid is at rest Ui = 0. Then, the 
component T xy is read off from (4.18), 


T xy = - Ph xy + P jicu?? + ^ (J 2 + q 2 ) k + ^iuj(u 2 - q 2 )p 

(uj 2 -q 2 )d- J(w 2 - g 2 ) 2 £j h xy . 
So, the retarded Green’s function in the scalar channel is 

G X ft' xy = — P + P jin;?? + ^ (uj 2 + q 2 ) k + |iw(w 2 — q 2 )p 

~\iu (^ 2 - Q 2 ) 0 - \(“ 2 - g 2 ) 2 £j . 


(E.5) 


(E.6) 


Shear channel It is sufficient to consider the metric perturbation ho x (v,z), which stirs 
the fluid to be of velocity u x (v,z). From the Navier-Stokes equations (E.2,E.4), 




iuiq 2 (7uj 2 — 3 q 2 ) 
192 



ho x 

—iui — q 2 t\ 


(E.7) 


Therefore, the component T 0x is 


T 0x = (e + P) u°u x + Pg 0x = (e + P) (- h 0x + u x ) + Ph 0x (E.8) 


which results in the retarded Green’s function in the shear channel 


G 0x ’ 0x = {e + P) 


8 q 2 i) - 2iuq 2 n - q 2 (q 2 - 2 cj 2 ) p - 4uj 2 q 2 6 - 2 iujq 2 (q 2 - uj 2 ) £ 
32 (— iuj + q 2 r]/ 4) 


-e (E.9) 


where the pole — iuj + q 2 g/ 4 = 0 is the shear dispersion relation. 


Sound channel The metric perturbation which generates the sound mode is ho z (v,z). 
The conservation laws (E.2,E.4) force a relation between the fluid velocity u z (v, z ) and the 
metric perturbation 


u z = <q + 


-iujq 


1 

64 


iujq 2 (7 uj 2 — 3 q 2 ) — 3q 2 t 3 


1 


16 x 36 
1 


iuj(21uj 2 + 31g 2 ) + (t-j — 2q 2 t.4) 
ho z • 


(E.10) 


q 2 — 3cu 2 — iujq 2 p + iuoq^C,/ 2 
Thus, the component T 0z is computed as 

T 0z = ( e + P) u°u z + Pg 0z = (e + P) (-h 0z + u z ) + Ph 0z , 
from which we read off the retarded Green’s function in the sound channel 


G 


0z,0z 

R 


= (e + P) 


IV* 


q 2 — 3c J 2 — iuq 2 r] + iu)q 4 (/2 


— £. 


(E.ll) 


(E.12) 
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The numerator IV* is 


at 2 -2 I ■ 4/- 1 2 2 1 ■ 3 2 I 1 • 3 2 n ^ 2 2 / 2 o 2'\ a 

IV* = q — zwg r/ + zccg C — 9 K ~ -«w g p + -iu; <7 v ~ ^ 9 w 9 (.9 — 3u; J 4 . 

From (E.12), we recover the sound wave dispersion relation g 2 — 3a ; 2 — iujq 2 r] + iujq 4: C > /2 = 0. 
Up to normalization convention and the 0-terms, our expressions for two-point correlation 
functions are the same as those of [3]. We also reproduce the plots [3] (also [51]) for the 
thermal correlators using our numerical results for the viscosities and GSFs. 

Choosing the spatial momentum along ^-direction, i.e., k tl = (—u), 0,0, < 7 ), the compo¬ 
nents of the correlators are conveniently written as 


Q&y,xy 


^ ,0a;,Ox 
G R 


£,00,00 

R 


q0z,0z 

R 


-G 3 (u,q), 


1 q 


=2 


2 a ) 2 - q 2 
2 


3 (a > 2 — q 2 ) 2 

2 Cd 2 q 2 

3 (w 2 - g 2 ) 2 


(a;,g), 

G 2 (a), < 7 ), 
C 2 (a), < 7 ), 


(E.13) 


where the bar is used to emphasize that the momenta here are dimensional. At zero 
temperature, the correlators in three symmetry channels coincide G\ = G 2 = G 3 = G s - 
For AT = 4 super-Yang-Mills theory, their expressions can be found in [51, 52] 

N 2 k 4 

G s = ( ln l fe2 | - i7r ® ( _fc2 ) si g n w) 1 (E.14) 

where the color N c is related to Gjy by 1/Gjy = 2 IV 2 / 7 t. We use our numerical results 
for viscosities and GSFs to construct the thermal correlators and find that they start to 
approach their zero temperature results (E.14) near a; q . 


F Perturbative solutions for V)’s and Tj’s 

In this appendix, we present perturbative solutions for Vi’s and Tj’s. Recall the formal 
expansion for the functions Vi’s and Tj’s in (4.23), 


Vi(u,qi,r ) = ^2 X n v} n \uj,qi,r), T(o;,%,r) = ^ \ n T$ n) (u, q if r). (F.l) 


H/ 


n =0 


n=0 


In what follows, we collect the results according to different sectors. 
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I: [Vi, T\. V 2 , T 2 } 


F 1 (0) = 0, V} 1 ’ = -icor 3 , 


(i) _ 


T (°) _ I 

1 “ 4 


In 


(1+r2 > ( 1 1+r)2 -2a r ctan( r ) + . 


r—>oo 1 1 




-~^ + o 4 


/*r />oo 

y 2 (0) = -y x 3 dxJ 


1 ^(y) 


yd 


3y 


dy - 


y 4t > ^ : 

i— 5-00 1 o , 

-)• --r z + O 



If = —zcj 


iuidx 


r X° - X Ji 
i'jj ( ln2\ / 1 \ 

* _ Ir 4 ( 1 2~~ ) ^ ^ ( r® ) ’ 


dy [2y 3 d y Tl°\y) -y + Zy 2 T[°\y) 


r—>oo 


= 


7^(0) _ 

1 '2 — 


A2, ~L x3ix i * [t^‘%) + ^ o (i), 

_2v^ M + nvr (!< ) + 2 (0)fa) 


l r x b -x J 1 


r—>oo 




y 

l 

r 

1—»oo 


dy 


I/ 2 (1) = -J™ x3 dxj^ dy^dyT^iij) O 


t[ 2) = 


dx r 

X — X 5 Ji 


dy 


2iuy 3 d y Ti 1 \y) + 3ia;y 2 T 1 (1) (y) + 


^ (2) (y) ^i (2) (d) 


rr> °°> 19 ^ 4 {6y 2 + u? [67t - 7r 2 + 12 (2 - 3In 2 hr 



(F.2) 


II: {V 3 , T 3 , V 4 , T 4 j 


vf = v 3 (1) = 0 , 


rp(0) _ 

3 “ 4 


(1 + r 2 )(l + r) 2 

In -- 4 -2 arctan(r) + 7r 


+ +0 (^5 


(o) f r , , /°° rXT. 0) (y , 1 r ^oo 1 2 ^ /1 

F; Uj = / x 3 dx / y 3 ^ ; dy + - — r 2 + 0 - 

4 1 J x 3 y y 8 12 1 


'i 


Tl(0) _ _ 

J 4 — 


dx 


x J — x Ji 


3 y 8 

?h>i + nh%) + 2 ,o, ■ 

y 2 y 3 


dy 


^ > 6 ^ _ 4^ (1 + 7r + 21n2) + °(^)’ 


7“>(1) _ 

J 3 — 


iudx 


x J — x 


dy 


2 y 3 d y T^\y) + 3 y 2 T 3 (0) (y) 


r ^oo iw io; . , . „/l\ 

—-t (—1 + In 2) + CM -7 ] . 

4r 2 8?’ 4 V r 5 / 


(F.3) 
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Ill: {V 5 , T 5 } 


Vi a) = -\r\ V 5 m =0, 


,(0) 1 , 1 + r 2 r-»oo 1 

' = - ln —2-> — 


T, w = 7 ^ 
4 


ri 1} = - 


r z 

00 icudx 


± 1 n ( 1 

* A„2 q„,4 ^ \ 


x 5 — X 


4r 2 8?’ 4 " T " ^ \r^) 

dy [2y 3 d y T^\y) + 3y 2 4 0) (y) 


IV: {T 9 } 


t 9 (0) = 0, 


^(l) 1 • 

T 9 = "I ZW 


™>^ ( 4 + ,-2h 2) + ofiV 

12 r 3 32r 4 \r 5 / 


(1 + r 2 )(l + r) 2 

ln --^-2 arctan(r) + 7r 


r ^oo zcu iu ( 1 \ 

+ I^ 4 +0 U)’ 


7^(2) _ 

^ 9 — 


= - [ t' Z _ [ d y ^y 3 d y T^\y)+3iuy 2 T^\y)+q 2 i 
Jr % % J 1 


r->oo <T +CJ 


V,VI,VII: {Vfj, T 6 , V 7 , T 7 , Vs , T s } 
V 6 (0) = 0, V 7 (0) = 0, 


^2“ + ^4 [g 2 + ( ln2 -!)w 2 ] +°(^) • 


T fi (0) = -T 7 (0) = * ln ^ 


K (0) = 


4 r 

/*oo 


_L , n(- 

4 r 2 g r 4 ^ ^ \^ r 5 I ’ 


r^f*vfM + i-oo.1 _a 

Ji Jr 3 y 


12 


*7^ + 0 

18 V r 


T“>(0) _ 

4 8 — 


/•°° da r 
Jr x 5 - x J l 


dy 


2 V 8 (0) (y) 2d y V 8 (0) (y) 

„,2 “ l " 




+fixffe) 


r^oo ) 13 - 12 In 2 + log r + ^ / logr \ 


oqc^4 


(F.4) 


(F.5) 


(F.6) 


The power expansion of U’s can be directly read off from these integrals. In principle, one 
can extend the above perturbative analysis to arbitrary order in derivative expansion. Then 
one can write down recursive relations for transport coefficients among different orders in 
the derivative expansion, as done in [2], 
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